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Abstract 

We consider the relative stability of parallel and perpendicular lamellar layers on corrugated 
surfaces. The model can be applied to smectic phases of liquid crystals, to lamellar phases of short- 
chain amphiphiles and to lamellar phases of long-chain block copolymers. The corrugated surface 
is modelled by having a single g-mode lateral corrugation of a certain height. The lamellae de- 
form close to the surface as a result of chemical interaction with it. The competition between the 
energetic cost of elastic deformations and the gain in surface energy determines whether parallel 
or perpendicular lamellar orientation (with respect to the surface) is preferred. Our main results 
are summarized in two phase diagrams, each exhibiting a transition line from the parallel to per- 
pendicular orientations. The phase diagrams depend on the three system parameters: the lamellar 
natural periodicity, and the periodicity and amplitude of surface corrugations. For a fixed lamel- 
lar periodicity (or polymer chain length), the parallel orientation is preferred as the amplitude of 
surface corrugation decreases and/or its periodicity increases. Namely, for surfaces having small cor- 
rugations centered at long wavelengths. For a fixed corrugation periodicity, the parallel orientation 
is preferred for small corrugation amplitude and/or large lamellae periodicity. Our results are in 
agreement with recent experimental results carried out on thin block copolymer films of PS-PMMA 
(polystyrene-polymethylmethacrylate) in the lamellar phase, and in contact with several corrugated 
surfaces. 
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1 Introduction 



Numerous ways to control the orientation of ordered meso-phases have been extensively studied in recent 
years. As an example, we mention that lamellar block copolymers (BCP) confined between two flat and 
parallel surfaces have been shown to orient parallel or perpendicular to the surfaces, depending on the 
surface separation Furthermore, chemically patterned surface can induce perpendicular or even 

tilted lamellae, if the chemical interaction is strong enough 00]. Applying shear is yet another effective 
method in producing large well- aligned samples but this method is difficult to implement in thin 
films. Another related situation can be found in liquid-crystals, where surface anchoring determines the 
direction of adjacent molecules and affects the bulk orientation and possible defects [7]. This type of 
surface effect has been studied extensively in relation to the Freedericksz transition and in twisted nematic 
liquid-crystal displays. Electric field can be very useful in aligning samples in which a large dielectric or 
conductivity contrast exists between the components |51 1§1 1101 1111 1121 IT31 1141 IT51 116j . However, the use 
of electric field requires rather sophisticated experimental setups in order to avoid adverse effects of ion 
accumulation at electrodes, short-circuits due to dust particles, etc. |16) . 

Most of the studies mentioned above are restricted to films in contact with flat and smooth solid 
surfaces. Little attention has been paid to the role of surface roughness on film morphology and ori- 
entation, although it is clear that non-flat surface topography gives rise to defects (e.g. vacancies in 
the lamellar ordering) and affects the lamellae orientation. Quite recently, experiments carried out by 
Hashimoto and co-workers on lamellar BCP films 17 addressed the question of how surface roughness 
affects BCP film orientation. In their study it was shown that the degree of surface roughness controls 
the lamellar orientation, leading to situations where the BCP films orient themselves parallel or perpen- 
dicular to the surface. These experiments serve as a starting point of our theoretical investigation, where 
we restrict ourselves to lamellar multilayer systems in contact with one corrugated surface. The model 
is expressed in terms of the elastic energy of lamellar or smectic systems, and is described by only two 
elastic constants. Hence, although some system-specific details are missing, the results are not restricted 
to BCP films, but are more general and equally apply to a broad class of systems ranging from smectic 
liquid-crystals to lyotropic (oil/water/amphiphile) systems. The main difference between these systems 
is in the values of the system parameters: the lamellar periodicity, strength of surface interaction and 
elastic constants. 

The elastic deformation energy of lamellar layers is studied separately for parallel and perpendicular 
orientations. It is shown that depending on the surface corrugation amplitude and periodicity, phase 
transitions can occur between the two orientations. Since real surfaces are never ideally flat, under- 
standing and characterizing surface roughness can be of great importance in controlling orientation of 
lamellar phases. This orientation mechanism is complementary to the mechanisms mentioned above of 
electric fields, shear and chemical surface patterning. In Sections [2 and [3| we calculate the deviations of 
lamellae from their corresponding flat perpendicular and parallel states. The free energy of the parallel 
and perpendicular states are compared in Sec. and a general discussion and comparison with previous 
works follows in Sec. 

2 Perpendicular Layers on Corrugated Surface 

Consider a lamellar system confined by one topographically corrugated surface, as depicted in Fig. 1. 
Instead of considering a rough surface having a random and quenched topography, we assume that the 
surface is characterized by a typical corrugation with a single amplitude and wavelength. The difference 
between a true rough surface and a corrugated one is not expected to be very significant as long as 
the rough surface amplitude and wavelengths do not vary much about their average values. A further 
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simplification is that the surface height h (measured along the z-direction) is taken to depend only on 
the lateral x direction, while it is translational invariant in the y direction: 

h(x) = Rcos(q s x) + ho , (1) 

where the average h = (h(x)) is taken hereafter to be zero, ho — 0, R is the corrugation amplitude and 
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Figure 1 : Schematic illustration of the rough confining surface. 

d s = 2-K/q s is the surface lateral periodicity. 

So far we described the quenched corrugated surface geometry. Next we consider the energetics and 
structure of a lamellar phase in contact with such a surface. We start by examining the order-parameter 
corresponding to lamellae which are perpendicular to the average surface position: 

4>±(r) = 0o cos(q x + qou{x, z)) , (2) 

where <po is the amplitude of lamellar concentration variations, d = 27r/g is the bulk lamellar spacing, 
and u(x, z) is a slowly-varying phase describing the deviation from flat lamellae, perpendicular to the 
surface and given by 4>o cos(qox). 

The total free-energy F of the lamellar stack in contact with one surface can be written as a sum of 
two terms, F = Fi, + F s , where F s is the surface energy and Fb is the bulk lamellar contribution. In a 
lamellar system there are different energy costs associated with bending and compression of the layers. 
It can be shown that for a slowly varying phase u(x, z\ an expansion of the free energy F\> to quadratic 
order in u and its spatial derivatives, can be written in complete analogy with smectic phases of liquid 
crystals as: 
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K {u zz f + B {u x f dV . (3) 



In the above equation the integral is over the entire volume, B is the compression modulus and K is 
the bending modulus. In BCP systems the model holds well in intermediate and strong segregations, 
where the two copolymer blocks are well separated in different domains and the chains are highly 
stretched. The elastic moudlii are K ~ do7AB and B ~ 7ab/<^0: where 7ab is the interfacial tension 
between the A and B polymer blocks, and do is the lamellar periodicity [T^]. For typical di-BCP such as 
polystyrene/polymethylmethacrylate (PS/PMMA), where 7ab — 2 mN/m and do ~ 50 nm, we estimate 
their values to be K ~ 10 -11 J/m and B ~ 4 • 10 5 J/m 3 . For liquid crystals, K ~ 2 ■ 10 -11 J/m and 
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B ~ 10 7 J/m 3 whereas for lyotropic (water/surfactant /oil) systems, K is about 1 fcsT/<io — 3T0 -30 
J/m and £? ~ 1 J/m 3 [20]. Clearly these lyotropic phases are very 'soft' and their elastic modulii have 
small values. 

The penetration length arising from Eq. (jSJ is 

A = {K/B) 1 ' 2 . (4) 

In lamellar di-BCP this length is proportional to the lamellar spacing, A ~ do — 50 nm. In lyotropic 
systems this distance is very small because of the reduced rigidity. 

The second term in the free energy is the surface contribution. Similar to previous works |4*l ll2lH5| . 
we assume a short-range surface field coupled linearly with the lamellar order-parameter at the surface. 



F s = J acj> dS . (5) 

The parameter a is the surface field and the integral is taken over the entire corrugated surface. The 
value of a is taken to be a constant throughout the surface, describing a corrugated surface which is 
chemically homogeneous. A positive (negative) a favors adsorption of negative (positive) (f> at the surface. 
In the language of an A/B di-BCP, this means that B (or A) monomers are preferentially adsorbed on 
the surface for a > (cr < 0). Assuming small distortions in the order parameter close to the surface, 
F s can be expanded to first order in u 

F s ~ 0-00 J [cos(q x) - Qou sm(q x)} dS (6) 

The entire profile 4>{x, z) can now be calculated by using a variational principle on the bulk free- 
energy, Eq. ©. In terms of the phase u(x, z), the resulting Euler-Lagrange differential equation away 
from the corrugated surface is, 

X 2 u zzzz — u xx = (7) 

This variation has to be complemented by a set of rather complex boundary conditions. They are 
obtained by taking the variation of the full F}, + F s on the corrugated surface, z — h(x), defined in 
Eq. CJ, 



df b 

du zl 



df s ( df b d df b 



du V dVu dz du z 



(8) 
(9) 



The unit vector normal to the surface is defined as: n = —(q s Rsinq s x, 1)/ ^(q s R sin q s x) 2 + 1, and fb 
and f s are the integrand of the volume and surface integrals, Eqs. and ©, respectively. 

In order to proceed we need to make some further simplifications. Solving the partial differential 
equation, Eq. (J7J, with the complex boundary condition, Eqs. is an extremely difficult task. A 

further simplification is to assume that u(x, z) is given by the single bulk mode <7o in the x direction, 
while it contains a sum over all possible Fourier modes in the z direction: 

u(x, z) = e iqoX A kt lkz + c.c. (10) 

k 

This simple sinusoidal form of u along the surface x direction may not fully account for the incommensu- 
rability that exists between the lamellar and surface periodicities. Defects and vacancies in the lamellar 
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ordering along the surface are not included either. Equations (Q and I|1U|) lead to a selection of specific 
k modes depending on the bulk mode go an d the penetration length A: 

k A + \- 2 ql=Q (If) 

This equation has two different roots, k± (and their two complex conjugates) 

k ± = ±(^y /2 .e±™/* (12) 

where exp(i7r/4) = (f + i)/\/2 is the 4th root of unity in the complex plane. The first boundary 
condition, Eq. JSJl, implies that Ak + — A^_ . Throughout this paper we assume q s R < qoR <C f , whence 
it is valid to approximate y/ (q s Rsmq s x) 2 + 1 ~ f and exp(ikR cos(q s x)) ~ f. The second boundary 
condition [Eq.^J] then gives 

-a<j)oqosmqox-iBq s Rsinq s x-J2loA k e t{ ' !oX+kz) - iK -^k 3 A k e t(q ° x+kz ^ = (13) 

k k 

The ^-dependent term in the exponentials is neglected as can be justified for q$R <C 1. The second term 
(proportional to B) can also be neglected for q s R <§; f , leading to the final expression of the lamellar 
order parameter: 

u = A e t( - qox+k + z ' ) + A Q e li - qax+k -^ + c.c. (14) 
~ 2K{k\ + kl) 

where A qo is denoted Aq in order to simplify the notation. The full order-parameter expression is obtained 
by substituting u from above in Eq. @ . The decay length of lamellar undulations is proportional to the 
lamellar spacing, l/k± ~ (A/go) 1 / 2 ~ do- 

The perpendicular lamellar stack on a rough surface is shown in Fig. 2. In Fig. 2(a) the lamellar 
periodicity is half of the surface one. The surface is attractive to the component marked in dark shades, 
causing them to expand in its vicinity. Light lamellae are in turn contracted close to the surface; this 
behavior is seen in part 2(c) too, where the surface periodicity matches the bulk lamellar one. Clearly, 
the curvature of the lamellae adjusts to the surface one in order to achieve the best compromise between 
elastic deformation and surface coverage. 



3 Parallel Layers on Corrugated Surface 

We now turn to describe a parallel lamellar stack in contact with the same sinusoidally corrugated 
surface. Our derivation is related to previous treatments [7||^|221j but the boundary conditions are 
handled differently by introducing the same type of surface term in the free energy F s , as in the previous 
section. The lamellae can be described along the same lines as for the perpendicular case but keeping 
in mind the important points where the stack orientation affects the free energy: 

</>l|(r) = -(f> cos(q z + qoutx, z)) (15) 

The perfect parallel layers, 0ii(r) = —<j) q cos q$z are recovered far from the surface (where u = 0). The 
bulk lamellar phase free-energy is obtained simply by interchanging the roles of x and z axes in the free 
energy of Eq. (3): x <-> z 

dV (16) 



F h = - 



K {u xx f + B (u z f 
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Figure 2: Perpendicular lamellar layers in contact with a sinusoidal surface, as obtained from Eqs. J2J) 
and JJ4|) . All length are scaled by do, the lamellar bulk periodicity, and the origin of the z axis is shifted 
by do/ 2 for clarity purposes, (a) The lamellar periodicity do is chosen to be half of the surface one, 
do = d s /2. Different grey shades correspond to different values of the order parameter cf> of Eq. J3), 
where dark and white correspond to negative and positive (f>, respectively, (b) Inter-material dividing 
surface corresponding to lines where 4> = in part (a), (c) The same as in (a), but the lamellar periodicity 
equals the surface one, d s — do- The parameters used are: B = 4, K = _B/(4gg) and a = \/BK/4:. 



As for the case of perpendicular lamellae, the surface energy, Eq. JSJ, can be expanded in small u as 
follows 



F s ~ (700 J [qousm{qoz) - cos{qoz)] dS (17) 

The governing equation for u is obtained from a variation principle applied to Eq. 116|) . SF^/Su = 

^?u xxxx — u zz — (18) 
with the boundary condition obtained from a variation of Ft, + F s on the z = h{x) corrugated surface. 



du 



( df d df 



o 



(19) 



\ <9Vit dx du xx 

Writing u as u(x,z) = Aoe*" 2 cos q s x and using Eq. I |18| l we find that a = iXq^. At this point it 
should be emphasized that when the surface periodicity is larger than the lamellar periodicity (q s < go); 
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the decay length 1/a ~ qo /q 2 for the parallel stack is much larger than the decay length ~ 1 /qo in 
the perpendicular case. This, in turn, means that for fixed R the elastic deformation gives preference to 
perpendicular ordering. 

As in the case of perpendicular lamellae, we use q s R < qoR <C 1 to approximate 
\J (q s Rsinq s x) 2 + 1 ~ 1 and obtain from the boundary condition, Eq. I|19|) . 

afoqlR + BXq 2 A = (20) 

yielding the order-parameter of parallel lamellae 

u(x,z) = A Q e taz cosq s x (21) 

° 5Ag? VBK\q s ) 

The parallel layering given by Eq. I|21() is plotted in Fig. 3. In 3(a) [as in Fig. 2(a)] the surface 
periodicity is twice larger than the lamellar one, and in 3(b) the two periodicities are equal. In the 
former case, dark regions (negative <fi) appear close to the wall, the lamellae are able to closely follow 
the surface contour and distortions are long-range. In the latter case the surface topography changes 
too quickly for the lamellae to follow and the lamellae lie almost perfectly flat. Distortions of the stack, 
in this case, can be seen only in the very close vicinity of the surface. 



4 Phase Diagram 

The elastic energy of the lamellar phase in the two orientations is obtained by substitution of Eq. (|14() 
in Eqs. 10 and JSJ (perpendicular lamellae) and Eq. (|21() in Eqs. I|16|l and JSJ) (parallel lamellae). 
Depending on the system parameters R, q s , qo and a, the minimum of the free energy is obtained for 
either one of the two orientations. 

Before we proceed in presenting the corresponding phase diagram, we note that our derivation 
is not valid over the entire i?, q s , qo and a parameter space. The assumption of small distortions, 
qou <C 1, together with Eq. (|21|l implies that Racftoq^qJ 2 <C s/BK. This can be rewritten as q s R 3> 
(qoR) 3 ^ 2 (cr(j)o/\ // B~K) 1 ^ 2 . Combining these two inequalities, we obtain that the limits of validity of our 
derivation are given by 

1 » q Q R > q s R » (<z i?) 3/2 (^L) (22) 
Cast in different terms this can be written as 

^ /2 ((?0 i?)i/ 2 «^ <1 (23 ) 

qo 



Before presenting the calculated phase diagram, it is of use to present a rough estimation of the 
parallel and perpendicular free-energies. Disregarding numerical prefactors, the bulk free energy of the 
perpendicular state is 



T ~ ^+^)-Im(b 
K q 



(24) 
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Figure 3: Parallel lamellar layers in contact with a corrugated surface, calculated from the free energy 
minimization of Sec. 3. (a ) The lamellar periodicity is half of the surface one do — d s /2, and the lamellae 
follow the surface topography, (b) The two periodicities are equal, d s — do, and the lamellae are unable 
to follow the surface height variations. Parameters here are the same as in Fig. 2. 

where S is the surface area and Aq is taken from Eq. (|14|l . At distances from the surface greater than 
1/Im(fc±), the distortion u is negligible, and the factor 1/Im(fc±) represents the effective volume to area 
ratio of the integration. 

The free energy of the parallel state is similarly estimated to be 



S uv " s > Im(a) 



-) (<Zoi?) 2 (25) 



K qo \q. 

Here Aq is taken from Eq. I|21(l and the effective volume to area ratio is 1/Im(a). 
The difference between the two free energies is roughly proportional to 

F\\ -fioc (O) (q Rf - 1 (26) 

Eq. (|26|l describes qualitatively the system behavior. If q Q is fixed while R and q s can vary, it directly 
follows that for qoR > q s /qo the perpendicular state is favored, while for q$R < q s /qo the parallel state 
is favored. Similar relations hold when q s is fixed and R and qo are allowed to vary, or when R is fixed 
but q s and qo can vary. 



8 



1 




Figure 4: Phase diagram in terms of q s and R with fixed surface field a and lamellar periodicity 2-k fqo- 
For fixed value of q s (scaled by qo ), a horizontal scan of increasing R (scaled by 1 /qo ) will increase the 
elastic penalty of parallel layers, and favors perpendicular ordering. On the other hand, a vertical scan 
of increasing q s (while keeping R constant) limits the deformations in the parallel state to the surface, 
and favors parallel ordering, a = 2\jBK and all other parameters are the same as in Figs. 2 and 3. 

The above simplified approach describes qualitatively the system behavior and may be used as a 
'rule of thumb'. Furthermore, the surface energies and the correct numerical prefactors can be taken 
into account as well. Figure 4 shows the phase diagram in the (q s ,R) plane, for fixed qo and o. Not 
the entire shown phase diagram is within the range of validity discussed above, Eq. (|23|) . For small 
values of q s , the perpendicular layering is favored, because of the long-ranged elastic strain pertaining 
in parallel lamellae as compared to perpendicular ones. As q s increases, the strain in the parallel state 
becomes more restricted to the vicinity of the surface, until, eventually, the parallel state becomes more 
stable. At this transition point, the energetic gain of having a commensurate layer close to the surface 
overcomes the loss of elastic energy deformation. 

Different conclusion can be drawn for small values of surface amplitude R, which generally induce a 
parallel state. Keeping q s fixed and gradually increasing R means that the elastic energy of deforming 
parallel lamellae increases, while the surface interaction stays constant. Therefore, at a certain threshold 
value of R there is a transition from parallel to perpendicular ordering. For larger q s values this critical 
R value increases as well. 

Figure 5 is a phase diagram in the (R, qo) plane, with q s kept fixed. For a given BCP chain length 
(fixed qo), increase of R will also increase the elastic energy of deforming parallel lamellae, and promotes 
perpendicular layering. On the other hand, keeping R fixed and decreasing go (so that it becomes 
comparable to q s ) decreases the range of parallel deformation; thus, favoring the parallel state. 

5 Discussion 

Lamellar stacks of either liquid crystals (smectics), short chains amphiphiles or long-chain BCP undergo 
deformation as they try to adjust to the presence of a rough (corrugated) surface. The amplitude and 
spatial extent of inplane and out-of-plane deformations are different and the preference for parallel or 
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Figure 5: Phase-diagram in the (i?, go) plane. The surface held a and surface periodicity 2n/q s are kept 
constant. Horizontal scans of increasing R (scaled by l/q s ) while keeping go hxed, increases the elastic 
penalty of parallel layers and favors perpendicular ordering. Same is true for vertical scans of increasing 
go (scaled by q s ) with hxed R. a = 0.2\/ BK and other parameters are the same as in Figs. 2 and 3. 

perpendicular orientations depends on specific bulk and surface parameters. Which morphology prevails 
depends on the surface periodicity d s — 2ix / 'q s , surface amplitude R, lamellar periodicity do = 2ir/qo, as 
well as the surface field strength a. 

For lamellae oriented perpendicular to the surface we use a model that indicates how the stack 
deformation propagates from the surface into the semi-infinite bulk. Parallel lamellae arc studied using 
a modification of a model used previously for smectics phase. In both cases small surface corrugation 
is assumed such that we can use the limits: q s R <C 1, qoR <C 1. Our study is complementary to a 
numerical study by Podariu and Chakrabarti who studied in detail extremely thin films (thickness 
comparable or smaller than the lamellar thickness) of a lamellar stack. 

The system behavior in terms of the parameters R and q s is given in the phase diagram of Fig. 4. 
For fixed surface periodicity 2ir/q s , increase in the surface amplitude R leads to preference of perpen- 
dicular lamellae. On the other hand, keeping R fixed and increasing q s leads to preference of parallel 
lamellae. This is a consequence of the diminished decay length of surface-induced undulations. When 
the undulations of parallel lamellae are restricted to the surface vicinity, the energy penalty of the elastic 
defect can be small. Thus, the phase transition between these two states is described by a line in the 
two-dimensional phase diagram. Quite generally, as the surface interaction parameter a is increased, 
this line moves towards the R axis in such a way that the state of parallel lamellae occupies a larger 
region in the phase diagram. 

The transition from parallel to perpendicular orientation as function of R and go is given in Fig. 5. 
Increase of R deforms the parallel layers and generally promotes perpendicular ordering. On the other 
hand, increase of go while keeping R and g s fixed, implies a reduced range of parallel deformations and 
yield a preference for parallel lamellae. 

Rough surface can be used to obtain morphologies that are usually controlled by chemical means: not 
only the transition between parallel and perpendicular lamellae, but also tilted lamellar morphologies. 
Indeed, in the limit of very small lamellar periodicity (large go) and strong interfacial interactions a, 
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lamellae will appear locally perpendicular to the surface, and therefore tilt when the surface is not 
horizontal. Similar tilted lamellae have been predicted for a BCP melt confined by chemically patterned 
surfaces [30], as the system tries to match the lamellae with the stripes. 

The present work is motivated by and is of direct relevance to recent experiments of Hashimoto and 
co-workers on PS/PMMA symmetric di-BCP on rough surfaces |17| . In the experiments, perpendicular 
orientation was observed for system parameters estimated to be q s ~ 0.04 nm -1 , qo — 0-33 nm -1 
and R ~ 7.5 nm. This set of data is in complete accord with the scaling formula, Eq. (|26[) . because 
q^R 2 /is — 1 is large and positive. This also can be confirmed by the phase diagrams presented in Figs. 4 
and 5. 

For different surfaces and BCP films reported in Ref. J7j, q s ~ 0.018 nm -1 , qo ~ 0.19 nm -1 and 
R ~ 2.7 nm, and the film orients itself in parallel layers. The scaling formula, Eq. (|26|l yields here a 
small positive number, implying weak preference for perpendicular layers. The phase diagram also show 
marginal behavior, possibly preferential for perpendicular layers. This discrepancy can have several 
(yet unknown) origins. First, the exact value of the surface interaction parameter a is not known from 
experiments, and it may be different than the values chosen by us in Figs. 4 and 5. Furthermore, we 
employed several approximations in our calculations ignoring, for example, possible vacancies in the 
perpendicular lamellae and the temperature dependence of the elastic modulii. We also assumed that 
the system is semi-infinite in the z-direction, while in experiments the BCP film thickness is finite. It may 
be that the BCP/air free surface can induce islands or other types of surface-induced defects that alter 
the simple picture employed by us. The free-energy model holds for BCPs in intermediate and strong 
segregations, but is less adequate for weaker segregations, as is possibly the case in the experiments. 
Lastly, our theory can be improved by taking into account more realistic surface roughness, instead of a 
single surface q-mode. 

We hope that additional and detailed experiments of lamellar systems in contact with well charac- 
terized rough surfaces will shed more light on this problem, and will further motivate theoretical studies 
of these intriguing systems. 
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